Insights on the cuprate high energy anomaly observed in ARPES 
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■Abstract 



^vq |Recently, angle-resolved photoemission spectroscopy has been used to highlight an anomalously large band renormaliza- 
^ .tion at high binding energies in cuprate superconductors: the high energy "waterfall" or high energy anomaly (HEA). 
Q^The anomaly is present for both hole- and electron-doped cuprates as well as the half-filled parent insulators with differ- 
<^ ent energy scales arising on either side of the phase diagram. While photoemission matrix elements clearly play a role 
in changing the aesthetic appearance of the band dispersion, i.e. creating a "waterfall" -like appearance, they provide 
an inadequate description for the physics that underlies the strong band renormalization giving rise to the HEA. Model 
calculations of the single-band Hubbard Hamiltonian showcase the role played by correlations in the formation of the 
HEA and uncover significant differences in the HEA energy scale for hole- and electron-doped cuprates. In addition, this 
OJ approach properly captures the transfer of spectral weight accompanying doping in a correlated material and provides a 
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unifying description of the HEA across both sides of the cuprate phase diagram. We find that the anomaly demarcates 
a transition, or cross-over, from a quasiparticle band at low binding energies near the Fermi level to valence bands at 
higher binding energy, assumed to be of strong oxygen character. 

'Keywords: ARPES, quantum Monte Carlo, Hubbard, strong correlations 
MCS: 74.72.-h, 79.60.-i, 74.25.Jb, 71.10.Fd 
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Advancements in angle-resolved photoemission spec- 
troscopy, an important probe of electronic structure, jl| 
have impacted significantly the study of strongly corre- 
lated materials. High resolution experiments, at binding 
energies up to 1 eV and higher, made possible by these 
advances, have revealed the presence of a "waterfall" -like 
structure with a characteristic kink at intermediate bind- 
ing energies - the high energy anomaly (HEA) - in the dis- 
persion of high Tc superconductors. j2l4l2| Found univer- 
sally in hole-doped compounds, the characteristic HEA en- 
ergy scale is ~ 300 meV with similar a dispersion anomaly 
"observed in the half-filled parent insulators. [l3| In con- 
.trast to earlier low energy studies focusing on dispersion 
"kinks under ~ 100 meV, [14] interpreted as signatures of 
coupling to low energy bosonic modes, [l5l-[l7| the extrap- 
olated band bottom has a value larger than that obtained 
from band structure calculations [2, 3] and the energy scale 
associated with the anomaly would tend to rule-out cou- 
pling to similar bosonic modes as the origin of the HEA. 

More recently, anomalies have been found in electron- 
doped compounds at approximately twice the energy scale. 

Taken together with results from hole-doped and 
half-filled parent materials, these findings raise questions 
about the origin or mechanism of the HEA, given the ap- 
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parent dichotomy in energy scales between electron- and 
hole-doped compounds, and what role, if any, many body 
effects may play. Focusing primarily on the HEA in hole- 
doped compounds, a number of theories have been ad- 
vanced including spin-charge separation, Q in-gap band- 
tails, [181 spin polarons, [19| couplingto spin fiuctuations, 
[6, 20, 21] phonons, Jj or plasmons, j22| strong correlation 
or "Mott" physics, [2|, |23I426| and even extrinsic^ffects as- 
sociated with photoemission matrix elements. [7] 

Photoemission matrix elements clearly play a role in 
modifying the appearance of the HEA, complicating the 
analysis. The kink- or "waterfall" -like structure found in 
the first Brillouin zone (BZ) in experiments performed us- 
ing synchrotron sources appears instead in newer exper- 
iments performed using low photon energy, laser sources 
as a band "break-up" or cross-over between a shallow, dis- 
persing band crossing the Fermi level and a higher binding 
energy feature near the F-point. [8] Matrix element effects 
also infiuence the results of investigations in higher BZs 
which find a shallow band with a characteristic "Y" ap- 
pearance near the zone center, rather than a true "water- 
fall". [7] However, the shallow, dispersing band near the 
Fermi level, together with the cross-over at the HEA en- 
ergy scale seen in both hole- and electron-doped cuprates, 

ill, J.2] originates from intrinsic band renormalization 
effects and not extrinsic mechanisms that merely serve to 
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change the appearance of features with changing experi- 
mental conditions. [13 Those theoretical scenarios based 
on weak coupling to high energy bosonic modes 0, [gI, [2qI - 
22I , 2?! have an appeal based on the kink-like appearance 
of the HEA, recalling earlier efforts aimed at explaining 
the origin of the low-energy kink in cuprates. l3- 17| 
However, coupHng to these modes, such as spin fluctua- 
tions, which would generally satisfy the energy scale for 
the HEA in hole-doped compounds, fails to account for 
the dichotomy in energy scales between hole- and electron- 
doped materials. However, spin fluctuations should play a 
prominent role in the renormalization mechanism forming 
the shallow, dispersing quasiparticle band for either hole- 
or electron-doped systems. 

To investigate the influence that many-body effects, 
i.e. strong correlations, may have on the origin of the 
HEA we investigate the single-particle spectral function of 
the single-band Hubbard model, building upon and adding 
to the inform ation obtained from much earlier investiga- 
Using quantum Monte Carlo [3l|,[32[ and the 
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tions. 

maximum entropy method (MEM) for analytic continua- 
tion, 33I, Q we study the spectral function for various 
values of electron filling. Our results indicate that the 
HEA can be connected to doping and the accompanying 
transfer of spectral weight into the lower or upper Hubbard 
band (LHB or UHB) of hole- or electron-doped systems, 
respectively. Doping leads to the formation of a quasi- 
particle band (QPB) at energies near the Fermi level, here 
set equal to zero energy, and the HEA represents a cross- 
over from this band to the LHB, playing the role of valence 
bands in the cuprates that would have substantial oxygen 
character. Correlations and spectral weight transfers in 
the model lead to a natural asymmetry in the HEA en- 
ergy scale between hole and electron doping, in agreement 
with experiment. The results also show qualitative sim- 
ilarities between the momentum dependence of the HEA 
from model calculations and that found experimentally. 

The single-band Hubbard Hamiltonian, an effective, 
low energy model for the cuprates derived from down- 
folding multi-orbital models explicitly incorporating cop- 
per and oxygen degrees of freedom, [sHI, [36| has the form 



where c] ^{ci^a) creates (annihilates) an electron with spin 



(J at site i, the site occupation n^^^ 



for each 



spin species equals or 1, {tij} is a set of tight-binding 
coefficients parameterizing the electron kinetic energy (in 
this study we restrict ourselves to nearest-neighbor, and 
next-nearest-neighbor, t' ^ terms), the chemical potential /i 
controls the electron filling, and the Hubbard repulsion U 
controls the strength of electron-electron correlations. The 
model is studied using determinant quantum Monte Carlo, 



31|, |32i] an auxiliary- field technique. This method yields 
the finite temperature, imaginary time propagator Gij{r) 
on a finite-size cluster with periodic boundary conditions. 




Figure 1: Bandstructure of the undoped single-band Hubbard model 
along high symmetry directions in the BZ. Model parameters: t' = 
-0.30t, = O.OOt, U = S.OOt. The UHB has a pronounced dispersion 
along the (tt, 0) — (tt, tt) and (0, 0) — (tt, tt) directions. The LHB, broad 
and centered at the F-point, has similar dispersing features, better 
separated from the bulk of the LHB especially along the (0, 0) — (tt, tt) 
direction. Note the well separated intensity near (7r/2, 7r/2) between 
— It and —2t and that at (tt, 0) near —2t. These dispersing features 
are precursors to the QPB that forms upon electron or hole doping. 
(Color online) 



Estimates for this propagator in imaginary time come 
from individual Markov chains of the Monte Carlo process. 
This imaginary time propagator must be Wick rotated, or 
analytically continued, to real frequencies to extract the 
spectral function. Performing a discrete Fourier trans- 
form and treating the ensemble {Gk(t)} obtained from 
individual Markov chains as a statistical sampling of the 
finite temperature propagator, the MEM, [33|, |3J based 
on Bayesian inference, is used to obtain the single-particle 
spectral function ACK^uj) on the appropriate discrete mo- 
mentum grid. 

The finite-size, square clusters used in this study have 
linear dimension N = 8 (64-site clusters). The correspond- 
ing momentum space grid {K} has a point spacing of 7r/4. 
The imaginary time interval, partitioned into L "slices" of 
size Ar = runs from to P = 1/T, the inverse tem- 
perature, and t serves as the energy unit of the problem. 
In this study, f3 = 3/t and Ar = 1/16 1. Once A(K,cj) 
has been obtained on the discrete momentum grid, the 
single-particle self-energy E(K,a;) can be extracted using 
Dyson's equation and the bare bandstructure correspond- 
ing to the tight-binding, model parameters. Assuming a 
weak momentum dependence to the self-energy, an inter- 
polation routine provides the value of E(k, a;) at an ar- 
bitrary point k in the BZ and Dyson's equation can be 
employed to compute A(k^uj). 

Figure [1] shows the bandstructure along high symme- 
try directions in the BZ for the undoped single-band Hub- 
bard model with parameters = — 0.30t, fi = O.OOt, and 
U = S.OOt. At the F-point and (7r,7r) one finds the bulk 
of the incoherent LHB and UHB, respectively. Above the 
Fermi level, the UHB has dispersing branches along the 
(tt, 0) — (7r,7r) and (0,0) — (7r,7r) directions with signifi- 
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Figure 2: Bandstucture of the hole-doped (~ 14%) single-band Hub- 
bard model along high symmetry directions in the BZ below the 
Fermi level. Each face represents a different symmetry direction ar- 
ranged in the proper geometric relationship to one another. The 
top face shows the Fermi surface with a portion removed, exposing 
the bandstructure along the BZ axis (the (0, 0) — (0, — tt) direction) 
and BZ diagonal (the (0,0) — (7r,7r) direction). Model parameters: 
t' = -0.30t, II = -2.50t, U = S.OOt. In this figure, the spectral func- 
tion A(k, ij) has been multiplied by the Fermi distribution function 
f(uj). The dispersing QPB crosses the Fermi level and remains shal- 
low, at low binding energies, giving way to the LHB, localized near 
the F-point, with a cross-over of weaker intensity at intermediate 
binding energies. (Color online) 




Figure 3: Momentum dependence of the hole-doped spectral function 
multiplied by the Fermi distribution function along selected cuts in 
the BZ. Model parameters: t' = -O.SOt, = -2.50t, U = S.OOt. 
Panels (a)-(d) show the intensity (falsecolor scale) along the cuts 
highlighted in the right panel. This sequence of cuts highlights the 
distinction between the QPB and the LHB and shows the evolution 
of the "waterfall" -like appearance of the HEA moving along the BZ 
axis away from the F-point. (Color online) 



cant spectral weight in the region near (7r,0). Below the 
Fermi level, the bulk of the LHB spectral weight is concen- 
trated near the F-point; however, there is less pronounced, 
but better separated, dispersing spectral weight along the 
(0, 0) — (tt, 0) and (0, 0) — (tt, tt) directions. 

There is a significant spectral weight in this dispersing 
feature near (7r/2,7r/2) between —It and —2t before los- 
ing intensity upon approaching the F-point and the bulk 
of the LHB. A similar energy scale can be assigned to the 
dispersing feature near (7r,0). The behavior of the spec- 
tral intensity along the (0,0) — (7r,7r) direction is qualita- 
tively similar to that observed in experiment [13] with a 
low binding energy feature near (7r/2, 7r/2) crossing-over to 
the valence band at the F-point with the transition marked 
by the appearance of a "waterfall" at intermediate binding 
energies. These dispersing features are precursors to the 
QPB that appears upon electron or hole doping. 

Figure [2] shows the bandstructure for the hole-doped 
single-band Hubbard model below the Fermi level along 
high symmetry directions in the BZ, arranged with the ap- 
propriate geometrical relationship. Here the spectral func- 
tion has been multiplied by the Fermi distribution function 
f{uj) = (exp(/3cj) + 1)"-^. The incoherent LHB, effectively 
localized near the F-point, has a weak tail of intensity ex- 
tending toward the points (tt, tt) and (0, — tt). The decrease 
in intensity within these tails approximately coincides with 
the location in momentum space identified with the HEA. 

Along the (0, 0) — (tt, tt) direction, the QPB disperses 
across the Fermi level at ~ (7r/2, 7r/2). Near (7r/4, 7r/4) the 
spectral intensity drops and the HEA appears as an appar- 
ent cross-over from the QPB to the LHB at an energy ~ 
-OM to -0.75t, the HEA energy scale. Along the (0, 0) - 



(0, — tt) direction, the QPB is nearly non-dispersive at an 
energy ~ — 0.5t. At approximately (0, — 7r/2) the QPB 
crosses-over to the LHB with a weak tail of intensity at 
low binding energy extending toward the F-point. While 
the spectral intensity decreases in the QPB approaching 
the F-point, there is significant coexistence of the LHB 
and QPB as a function of momentum along these high 
symmetry cuts. The spectral function has a "waterfall" - 
like appearance even without including the effect of pho- 
toemission matrix elements, which serve to enhance the 
appearance by further suppressing intensity nearing the 
F-point. [11] 

The relationship between the QPB, the LHB, and the 
appearance of the HEA can be further explored studying 
the detailed momentum dependence of the spectral func- 
tion. Figure [3] shows the spectral function along selected 
momentum space cuts. Cut (a) reproduces a part of the 
spectral function along the (0, 0) — (tt, tt) direction already 
encountered in Fig.[2j Cuts (b)-(d) show the spectral func- 
tion in momentum space, parallel to the (0, 0) — (tt, tt) 
direction, moving toward (0,7r). While the energy scale 
of the LHB remains relatively unchanged, its spectral in- 
tensity progressively decreases. The QPB, crossing the 
Fermi level, increasingly becomes easier to identify as a 
proper band with increased spectral weight approaching 
the BZ axis. The "waterfall" -like appearance of the HEA 
also evolves with changes to its momentum space position 
and progressive reduction in spectral intensity below the 
QPB. Taken as a whole, the evolution of the QPB and 
HEA qualitatively agree with the results of experiment on 
hole-doped compounds, [H, Q including the evolution of 
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Figure 4: Bandstucture of the electron-doped (~ 17%) single-band 
Hubbard model along high symmetry directions in the BZ below the 
Fermi level. The top face shows the Fermi surface with a portion 
removed. Model parameters: = -0.25t, = 2.20t, U = 8.00t. In 
this figure, the spectral function A(k, uu) has been multiplied by the 
Fermi distribution function f(uj). The dispersing QPB crosses the 
Fermi level and remains shallow, at low binding energies approach- 
ing the F-point with diminishing intensity. The LHB (not shown) 
at higher binding energies (~ — 7.5t) has a substantially reduced in- 
tensity due to spectral weight transfer into the QPB and UHB with 
electron doping. The HEA energy scale is ~ —l.Ot to — 1.5t, twice 
that found in the hole-doped system. (Color online) 
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Figure 5: Momentum dependence of the electron-doped spectral 
function multiplied by the Fermi distribution function along selected 
cuts in the BZ. Model parameters: t' = -0.25t, /i = 2.20t, U = S.OOt. 
Panels (a)-(e) show the intensity (falsecolor scale) along the cuts 
highlighted in the lower right panel. The lack of significant spectral 
intensity in the LHB suppresses the "waterfall" -like appearance of 
the HEA. (Color onhne) 



spectral intensity and changes in momentum space posi- 
tion and robustness of the "waterfan"-hke appearance as 
a function of momentum. 

Comparison of calculation results for a hole-doped sys- 
tem to those for an electron-doped system reveal the di- 
chotomy in HEA energy scales intrinsic to the single-band 
Hubbard model. Figure |4] shows the spectral function mul- 
tiplied by the Fermi distribution function for an electron- 
doped system with model parameters t' = — 0.25t, /i 
2.20/. and U = S.OOt. The LHB (not shown) has a dra- 
matically reduced intensity due to spectral weight trans- 
fers into the QPB and UHB upon electron doping. This 
affects the "waterfall" -like appearance of the HEA in the 
electron-doped calculation by reducing the spectral inten- 
sity at intermediate binding energies between the QPB 
and the LHB. Previous results demonstrate that modify- 
ing the spectral intensity through the inclusion of approx- 
imate matrix elements leads to enhanced "waterfall" -like 
character in the electron-doped dispersion. For real- 
world experiments, the proper valence band, composed 
of significant oxygen character, should remain robust un- 
der electron-doping, increasing the likelihood of a promi- 
nent "waterfall" -like appearance in the cross-over region 
at intermediate binding energies, in agreement with ex- 
periment, [ol-fiil 

The dispersive QPB shown in Fig. |4] dips farther be- 
low the Fermi level than does its hole-doped counterpart, 
yielding an HEA energy scale ~ —IM to —IM along the 
(0,0) — (7r,7r) direction. Along the (0,0) — (0, — tt) direc- 
tion, the spectral function appears to be quite similar to 
its hole-doped counterpart with an energy scale ~ — 0.5t 
near (0, — tt) with a weak dispersion in this region. To- 
gether with decreasing intensity, the QPB along this di- 



rection disperses downward to and energy ^ —IM at the 
F-point. The chief contrast between the QPB in hole- and 
electron-doped systems, other than the energy scale, lies 
in their origin. With hole doping, the chemical potential 
moves down into the LHB, or more precisely the dispersive 
features at lower binding energy, and the QPB originates 
from the dispersive precursors shown in Fig. [TJ whereas 
in electron-doped systems, the chemical potential shifts in 
the opposite direction and the QPB arises from the pre- 
cursors in the UHB. 

Figure [5] shows the momentum dependence of the QPB 
along selected cuts in the BZ. Like the hole-doped system, 
the spectral intensity near the BZ axis progressively in- 
creases moving from the F-point toward (7r,0) or (0,7r). 
While there is significantly less intensity below the QPB 
compared with the hole-doped system, one can also in- 
fer changes to the HEA momentum space position and its 
"waterfall" -like appearance from cuts (a)-(e) that would 
be qualitatively similar to the hole-doped system. The 
most noticeable contrast to the hole-doped system is the 
change in the HEA energy scale between cuts (a) and (e), a 
factor on the order of 2 or 3, seen directly in the cut along 
the BZ axis (cut (a) of Fig.[5|). As in the hole-doped case, 
these results qualitatively capture the behavior observed 
in experiments on electron-doped compounds. 
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The results presented in this study provide evidence 
for an HEA in the single-band Hubbard model. The find- 
ings appear to agree qualitatively (quantitatively with a 
proper choice for the energy scale t) with experiments on 
both hole- and electron-doped compounds as well as the 
half- filled parent insulators. In principle, the results would 
include effects associated with the coupling of electrons to 
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spin fluctuations, one of the theoretical scenarios for the 
origin of the HEA, at an energy scale proportional to the 
superexchange constant J that should have a similar value 
in hole- and electron-systems. However, the dichotomy in 
the HEA energy scale between hole- and electron-doped 
systems (a factor of ~ 2) would argue against solely this 
origin due to bosonic-mode coupling. Instead these results 
support the conclusion that strong correlations play a cen- 
tral role in the origin of the HEA and that the anomaly 
itself results from a simple cross-over between the shallow, 
dispersing QPB near the Fermi level and the incoherent 
LHB at higher binding energy. 
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